Are all modes created equal? An analysis of the WMAP 5- and 7-year data without 

inflationary prejudice 
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We submit recent claims of a semi-significant detection of primordial tensor perturbations in the 
WMAP data to a closer scrutiny. Our conclusion is in brief that no such mode is present at a 
detectable level once the analysis is done more carefully. These claims have their root in a brief 
debate in the late 1990s about the standard calculation of the scalar and tensor spectra in standard 
inflationary theory, where Grishchuk and collaborators claimed that their amplitudes should be 
roughly equal. We give a brief summary of the debate and our own reasons for why the standard 
calculation is correct. 

PACS numbers: 98.80.-k, 98.80.Cq, 95.85.Sz 



I. INTRODUCTION 



Cosmic inflation has become a cornerstone of what may 
be called the standard model of cosmology. Apart from 
its role in lessening the demands on the need for fine- 
tuned initial conditions for the observable universe, the 
major success of the inflationary paradigm has been to 
explain the existence and statistical properties of the per- 
turbations in the distributions of matter and radiation 
(see for a gentle introduction). Inflation, in the sim- 
plest versions realized by a self-interacting scalar field, 
predicts a nearly scale-invariant initial power spectrum 
of density perturbations. The subsequent evolution of 
these perturbations limits the amount of information we 
obtain from observations of the statistical distribution of 
matter in the present universe. This is a pity, since infla- 
tion most probably involves physics at energies exceed- 
ing those obtainable in particle accelerators by several 
orders of magnitude, and can therefore give important 
constraints on theories that go beyond the present Stan- 
dard Model of particle physics. 

Fortunately, another robust prediction of inflation is 
that it also blows up the quantum fluctuations in the 
spacetime metric to observable scales and produces a 
spectrum of gravitational waves 0. These remain un- 
affected by the subsequent cosmic evolution, and there- 
fore provide a window to the physics of the inflationary 
process itself. Unfortunately, the standard prediction of 
single-field inflation is that the amplitudes of these waves 
are much lower than those of the density perturbations, 
and are hence much harder to observe. This standard lore 
has, however, been challenged, at several times and in 
various papers, by Grishchuk (see e.g. 0,0]). His claim 
is that the standard inflationary treatment of perturba- 
tions is wrong, and that when the calculation is done cor- 
rectly what emerges are spectra of density (scalar) and 
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gravitational wave (tensor) perturbations with roughly 
the same amplitude. Furthermore, in a series of recent 
papers 0-0 he and collaborators claim to find weak ev- 
idence for a gravitational wave signal in the WMAP 5- 
and 7-year data 0, 0- If his claims are correct, a con- 
vincing detection of primordial gravitational waves may 
be very close, perhaps already with Planck data. They 
therefore merit closer scrutiny, which is what we set out 
to provide in the present paper. 

The structure of our paper is as follows: In section 
II we summarize and comment on the analyses of the 
WMAP 5- and 7-year data found in Zhao et al. 0, 0, 
before giving an outline of our strategy in section III. In 
section IV we present results from a full analysis of the 
WMAP 7-year data, taking all relevant parameters into 
account. This provides a background for judging the re- 
sults of the more restricted analyses of Zhao et al., which 
we carry out improved versions of in section V. Section 
VI contains our comments on the theoretical aspects of 
the disagreement on the relation between scalar and ten- 
sor perturbations. Finally, we summarize our findings 
and conclude in section VII. 



II. A SUMMARY OF THE ANALYSIS IN ZHAO 
ET AL. 

We will in the current section briefly describe the anal- 
yses performed in and , and the results and conclu- 
sions of those analyses, in order to motivate the analyses 
performed in the current paper. 

Both papers contain analyses of WMAP data and fore- 
casts for the Planck mission. Our concern in this paper 
is the data analysis part, and so we only present the 
methodology, results, and conclusions of the WMAP data 
analyses. 

The authors first derive the likelihood function, which 
depends on the true CVs as well as the observed data 
points and noise values. They use the WMAP data points 
and noise values as provided at the LAMBDA web page 
[30| . although they approximate the noise values to be 
uncorrelated between different multipoles. 
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Using the MCMC method, the authors use the 
WMAP5 (in @) and WMAP7 (in @) data to map the 
likelihood functions for three parameters: A s (the ampli- 
tude of the scalar perturbations), n s (the scalar spectral 
index), and R. The first two of these are the standard 
cosmological parameters, while R is the tensor-to-scalar 
quadrupole ratio, which the authors use instead of the 
standard parameter r, the tensor-to-scalar power spec- 
trum ratio at k = kj, = 0.002 Mpc" 1 . (These are roughly 
related by r ss 2R [5j.) The other cosmological parame- 
ters are fixed to their mean values from the WMAP team 



5-year Q and 7- year [l(| analyses. 

The authors only use the multipole ranges 2 < I < 
100 and 101 < £ < 220 in their analyses - the former 
for the actual analysis, and the latter in order to show 
that gravitational waves have little or no effect beyond 
I fa 100, and that it therefore is futile to use multipoles 
higher than this in the search for gravitational waves. 
Also, they only use the TT and TE power spectra, saying 
that the EE and BB modes are not very informative. 

For their 2 < I < 100 analyses, the authors find the 
following 3-dimensional ML values: 



R = 0.229, n s = l.Q86, A s = 1.920 x 10~ 9 (1) 

for the 5-year analysis and 

R= 0.264, n s = 1.111, A s = 1.832 x 10~ 9 (2) 

for the 7-year analysis. Further, when marginalising the distributions for each of the parameters, they find the 
following peak values and Minimum Credible Interval (MCI) (see [ll| for an explanation of this interval) limits: 

R = 0.266 ± 0.171, n s = 1.107±g;^, A s = (1.768±g j^) x 10~ 9 (3) 

for the 5-year analysis and 

R = 0.273+8:^, n s = 1.112™, A s = (1.765t°; 2 26 l) x 10" 9 (4) 

for the 7-year analysis, where all limits are the 68.3% confid ence intervals. 



These results contrast with those of the WMAP team, 
which finds no compelling evidence for gravitational 
waves 0. In order to explain this discrepancy, Zhao 
et al. suggest (as we have mentioned) that it is wrong 
to use higher multipoles than £ = 100 in the search for 
gravitational waves. In addition, they point to their 
101 < I < 220 analyses, where they find the follow- 
ing 68.3% MCI limits for n s : n s = 0.948±g;g|f (5-year- 
analysis) and n s — 0.969^Q;gg3 (7- year analysis). Com- 
paring with the £ — 2 — 100 analysis, these MCI limits do 
not overlap (5-year) or just barely overlaps (7- year) with 
each other. This suggests to the authors that the assump- 
tion of a £ (or fc)-independent spectral index is erroneous, 
and this assumption combined with the unwarranted use 
of multipoles larger than 100 is what leads the WMAP 
team to the wrong result. They also mention that the 
WMAP team uses inflationary predictions, in particular 
the consistency relation r — —n t /8 (where n t is the ten- 
sor spectral index), in their analyses, which makes the 
analysis biased from the outset. 



A. Preliminary comments 

We here make some short comments about the conclu- 
sions reached in the two papers by Zhao et al. 



First of all, we disagree with their claim that when 
looking for gravitational waves, one should only include 
the multipoles that are affected by them. As the authors 
themselves have pointed out, there are degeneracies be- 
tween several parameters at the low £ range: A s , n s , 
and R can be modified so that two different values of, 
say, R can produce almost exactly the same power spec- 
trum at large scales. Hence, it makes sense to include as 
much data as we can in the likelihood analysis: Small- 
scale data may perhaps constrain A s or n s , which in turn 
gives better constraints on R, since the degeneracies may 
then be broken to some degree. If, when including small- 
scale data, we get a small value for R, it simply means 
that there are very few models where a large R can work 
together with n s and A s (and possibly other parameters) 
to produce a power spectrum that fits both on large and 
small scales. 

Next, concerning the non-compatible values of n a for 
2 < £ < 100 and 101 < £ < 220: There does indeed seem 
to be a tension between these results and the assump- 
tion of a constant spectral index. However, a la signal 
is a hint, not a detection. Further, even if the assump- 
tion of a constant spectral index turns out to be wrong, 
it does not follow that this assumption is what leads to 
the discrepancy between their analyses and those of the 
WMAP team - the solution to the n s -tension could still 
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yield a negligible R, especially considering the degenera- 
cies mentioned above. More on this below. 

Finally, we agree with the authors that the analysis of 
the WMAP team may be biased by the explicit use of the 
consistency relation in their data analysis: This relation 
should appear automatically if the theory of single-field 
inflation is correct, and not be used as a constraint in the 
analysis. However, we also recognise the need to reduce 
the number of free parameters in the likelihood analysis, 
which the WMAP team cites as their reason for using 
the consistency relation. 



III. OUTLINE OF OUR STRATEGY 

Our aim in this paper is both to verify the analyses 
performed in [B| and [g] , but also to further test the con- 
clusions they draw from these analyses. To this end, we 
will utilise the MCMC method as Zhao et al. do, but our 
approach will differ somewhat from theirs. 

Most importantly, we will be using the CosmoMC soft- 
ware [l2j to perform the likelihood analysis. The advan- 
tages of using this rather than a self-made sampler should 
be obvious: It is well-tested, and bugs and errors are less 
likely to occur. The additional advantage is that this 
software is well integrated with the WMAP likelihood 
software, so that there is no need for e.g. making ap- 
proximations to the noise values as Zhao et al. do. All 
power spectrum modes can be incorporated without ex- 
tra effort as well. Also, it allows for simulation of the 
lensing effect on the power spectra, and we use this fea- 
ture throughout our analyses. 

To start off, in section ITVl we first do a full analysis of 
the WMAP7 data for 2 < £ < £ ma x, where £ max is the 
highest multipole for which there are data, and all power 
spectrum modes. In this analysis we use the constraint 
ris = rit — 1 instead of the consistency relation, following 
[J-Q where this constraint was imposed for theoretical 
reasons. This gives an indication of what the WMAP 
data would give when the consistency relation is not en- 
forced. 

We analyse the low I WMAP data in section IVl First, 
in section lV Al we repeat the analyses performed by Zhao 
et al. in order to see if the use of the exact noise values 
and CosmoMC software has any impact on the results. 
We choose to vary the optical depth to reionization, r, in 
addition to the three parameters mentioned above, since 
this parameter can be expected to share some degeneracy 
with the three others. We suggest a way to obtain more 
appropriate best-fit values at which to fix the parameters 
we do not vary, and repeat the analyses using these best- 
fit values. 

Taking our \ow-£ analysis a step further, in section IV Bl 
we aim to test the claim made by Zhao et al. that the 
assumption of a constant n s is what leads the WMAP 
team to overlook the gravitational wave contribution to 
the CMB. To do this, we introduce a 'step-like' spectral 
index with a jump at i = 100, and do analyses for the 



range £ = 2 - 220. 



IV. FULL ANALYSIS OF ALL WMAP7 DATA, 
WITH ALL RELEVANT PARAMETERS 

In this section, we do an initial full analysis of the 
WMAP 7-year data. We use, as already mentioned, the 
CosmoMC sampler to find the likelihood function for a 
7-dimensional parameter space, consisting of the parame- 
ters n h h 2 , rj D M^ 2 , t, n s , log(10 10 A s ), 6, and R. All these 
parameters are standard CosmoMC parameters (ft^h 2 is 
the physical baryon density, iluMh 2 is the physical den- 
sity of cold dark matter), except i?, with 9 being the 
ratio of the sound horizon to the angular diameter dis- 
tance, which is being used instead of the Hubble param- 
eter because it is less correlated with other parameters. 
We modify CosmoMC to use R instead of r. As also men- 
tioned before, we use the relation n t = n s — 1 instead of 
the inflationary consistency relation. Following Zhao et 
al. and the WMAP team, we use k p = 0.002 Mpc" 1 as 
the pivot wavenumber. 

This choice of parameters implies certain assumptions 
- namely, a flat Universe, dark energy with the equation 
of state w = — 1, massless neutrinos, and no running of 
the spectral index. 

We do not use a fixed number of samples; since we do 
8 chains simultaneously, we use the CosmoMC feature to 
check for convergence of confidence limits: By computing 
the variance between chains of 2.5 per cent of the distri- 
bution tails, the error must be smaller than 0.2 in units 
of the distribution's standard deviation. This means that 
we can be sure that the sample confidence limits are not 
too far away from the true confidence limits. 

We use the four power spectrum modes TT, TE, EE, 
and BB, and we do the analysis for the range 2 < I < 
£max, where £ m ax is 1200 for the TT mode, 800 for the 
TE mode, and 23 for the EE and BB modes. 



A. Results 

The results of the analysis described above are shown 
in table fl] Both full-dimensional ML values and one- 
dimensional peaks and MCI limits are shown. 

B. Discussion 

We have performed a full-scale likelihood analysis of 
the WMAP 7-year data, not using the consistency rela- 
tion, in accord with the complaints voiced by Zhao et al. 
Despite this, there is little reason to claim any detection 
of gravitational waves in the CMB, as the R = hypoth- 
esis is well within confidence bounds. This means that 
the only explanation left, as offered by Zhao et al., to 
explain the discrepancy between their results and those 
of the WMAP team must be the (erronously?) assumed 
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Parameter ML value One-dimensional peaks and MCI limits 



Y 68%t, 95%t 


n b h 2 


0.0227 


O 09^9 0.0239. 0.0248 
U.UZOZ q 0224, 0.0217 


"DM AI 


0.1107 


~ * 01130 11 8Q 
U.1UUU Q.1003, 0.0936 




1.040 


i n/n 1.044, 1.047 
1.038, 1.035 


r 


0.0880 


n n«Q/i °- 1051 . 0.1229 

u ' uotl * 0.0738, 0.0611 




0.970 


n q«9 1 007 > 1 034 

U.yOZ o.964, 0.947 


log(10 10 A 3 ) 


3.174 


o , n n 3.195, 3.249 
3 052, 2.958 


R 


0.012 


n n1 n 0.089, 0.174 
U.UlU o.ooo, 0.000 



TABLE I: ML points and marginalised one-dimensional peaks 
and MCI limits for a likelihood analysis for seven parameters 
using the WMAP 7-year data. MCI limits are given with both 
68% and 95% confidence. 



constancy of the spectral index for all £'s. We will test 
this explanation presently, but first we will repeat their 
analyses of the low I WMAP data. 



V. LOW-f ANALYSIS 

We will now focus on the large-scale WMAP data, as 
the central claims of Zhao et al. pertains to these scales. 
First, we repeat the work of Zhao et al. in section fV Al be- 
fore we move on to testing the authors' claim concerning 
the non-constancy of the spectral index in section fV Bl 

We will need to limit the number of free parameters 
to use in this section, as we will only be using multipoles 
up to I = 220, and this data range contains too little in- 
formation to be able to constrain all parameters at once. 
Which parameters we vary will be specified below. We 
will, however, never use the Hubble parameter, instead 
using 9. This then differs from the approach of Zhao et 
al. 

The parameters we do not vary must be given a certain 
value. As mentioned above, Zhao et al. have chosen to 
use the mean value of the marginalised one-dimensional 
distributions for each parameter from the WMAP stan- 
dard analysis. We find two problems with this approach: 
First, the WMAP standard analysis excluded gravita- 
tional waves, making the best-fit parameters obtained 
less valid when we want to include gravitational waves 
in our analyses. Second, we think that the mean values 
are not the best values to use. The one-dimensional peak 
values would be better, as these represent the most prob- 
able value of each parameter given the data we have (of 
course, for Gaussian distributions, the mean and peak 
values coincide) . In our opinion it is better to choose the 
parameter values corresponding to the peak of the full, 
multi-dimensional likelihood function: These taken to- 
gether are the closest to what one may call the Standard 
Model. 

In order to remedy these issues, we here do the fol- 
lowing: We do an initial analysis of most available data: 
WMAP 7-year data (2 < I < £ m ax), matter power spec- 



Parameter ML value One-dimensional peaks and MCI limits 



Y 68%t, 95%t 
A 68%J., 95%J. 


n b h 2 


0.0226 


n 099S 0-0233, 0.0238 
0.0223, 0.0218 


^DMh 2 


0.1177 


D 1 1 7fi 0.1203, 0.1229 
U.ll/O 0.1145, 0.1118 


e 


1.040 


1 040 1 043 ' 1045 

J..U1U L038, 1.036 


r 


0.0810 


n n«7R - 1028 . 0.1176 

u.uem 0.0749, 0.0625 


n s 


0.963 


n QfiO. 0-983, 0.996 

u.yoy 0.957, 0.945 


log(10 10 A a ) 


3.207 


n ons 3 - 243 > 3 281 
3 165i 3 12 3 


V 


0.015 


000 °' 076 ' 0168 

U.OUO 0.000, 0.000 



TABLE II: ML points and marginalised one-dimensional 
peaks and MCI limits for a likelihood analysis for seven pa- 
rameters, taking into account CMB, power spectrum, and su- 
pernova data, in addition to various prior constraints on some 
parameters, described in the text. MCI limits are given with 
both 68% and 95% confidence. 

trum data (SDSS, fourth data release [13); supernova 
observations (SDSS, ESSENCE [3, SNLS HST 
16] and various low-redshift supernova observations), 
Lyman-alpha data (SDSS), and priors on the age of the 
Universe, the Hubble parameter (from the HST), and fib 
(from BBN). We let vary the same parameters as in sec- 
tion II V| but use r instead of R. We do this as R is harder 
to implement for non-CMB data, since it is just the ratio 
of the tensor contribution to the CMB quadrupole to the 
scalar contribution, while r is more 'fundamental', as it 
describes the tensor/scalar ratio of the primordial power 
spectra. We do keep using the relation nt = n s — 1. We 
then take the full-dimensional ML values from this anal- 
ysis and use these as best-fit values for the rest of the 
analyses (though we will also use the best-fit values used 
by Zhao et al., see below). 

The results of this analysis are shown in table HU We 
will from here on carry out two versions of the analysis: 
One using the ML values found in table [Hi and one us- 
ing the best- fit values used by Zhao et al., in order to 
examine whether using different best-fit values makes a 
large difference in the results. When using the best-fit 
values from Zhao et al., we use the 5-year values for the 
5-year analysis, and 7- year values for the 7-year analysis. 
We will, as mentioned, not use the Hubble parameter, so 
in this we do not follow Zhao et al. The value for 8 will 
always be the ML value from table |Hj 



A. Analysis with CosmoMC 

In this section, we repeat the \ow-£ analyses in Q and 
with certain elaborations and modifications, which we 
here specify: 

As before, we use CosmoMC with the WMAP like- 
lihood software, and include all power spectrum modes 
in the analysis. We do analyses for both 5-year and 7- 
year WMAP data in order to compare with both papers 
by Zhao et al. We include r as a parameter to be var- 
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ied, though we do analyses where we vary only the three 
original parameters as well. We do all analyses over three 
ranges: 2 < I < 100, 101 < I < 220, and 2 < I < 220. 
(The £ = 101 — 220 contained too little information to 
constrain four parameters at once, so we did not include 
r in the analyses for this range.) We do an additional 
run where we fix R to be its full-dimensional ML value 
from the I = 2 — 100 analysis where r is not varied. This 
run is done for later comparison. 

1. Results 

The results of the above analyses are shown in tables 
IIIII f for the revised best- fit values) and llVl ffor the best- fit 
values used by Zhao et al.) We have assigned to each run 
a number for easier referral. 



2. Discussion 

From tables UTT1 and HVT we note the following features 
of our analyses: The values for R are consistently lower 
than those found in the analyses of Zhao et al., and the 
one-dimensional distributions of n s for runs 1 and 2 over- 
lap within a la interval. This seems to hold no matter 
which best-fit values we use (though with the best-fit val- 
ues used by Zhao et al. we find a slightly higher R than 
with the best-fit values derived here), which indicates 
that the difference in results from the analysis of Zhao 
et al. most likely is due to use of the WMAP likelihood 
software with exact noise values (The use of the EE and 
BB modes should have less of an impact due to the high 
noise in these modes). 

The inclusion of r actually seems to raise the values 
for R slightly, and we see what Zhao et al. also reported: 
R has higher values in the 7-year analyses than in the 
5-year ones (mostly). Also, in some cases, we do find 
lower limits for R. However, these are very small, and 
we only have a la lower limit, so it is hard to draw any 
conclusions from these results. 

The fact that the distributions for n s overlap within la 
between runs 1 and 2 weakens the claim that the spectral 
index really is dependent on I. Even though the overlap 
is slight, we must again be reminded that a la detection 
is very weak, statistically speaking. Nevertheless, in the 
next section we will try to put the claims about n s made 
by Zhao et al. to the test. 

B. Testing the hypothesis of an ^-dependent n s 

In this section, we wish to test, on a very basic level, 
whether introducing an ^-dependent spectral index gives 
better fit to data, and whether it may yield a high value 
for R over a larger multipole range than I — 2 — 100. 

In order to accomplish this, we implement a 'step-like' 
spectral index: The spectral index is allowed two values 



- rig 100 below I — 100 and n s >10 ° above this multipole. 
This is not implemented directly - rather, I = 100 is 
translated into fc-space through the relation [l7| (- ~ kd a , 
where 



is the angular diameter distance, and Q m — fluM + is 
the density parameter of dust. We choose to work in k- 
space both because it is easier to implement, and because 
it makes slightly more sense physically. However, we em- 
phasise that there is no physical justification for such a 
spectral index - we are just attempting to test a claim 
which will eventually need some physical justification if 
found to be probable. 

We still use the relation n t — n s — 1, so that both 
the scalar and tensor spectral indices attain this step-like 
behaviour. 

Another, more straightforward way to implement an 
^-dependence would be to use a simple running of the 
spectral index. However, we have chosen not to do this 
because this case is discussed in Zhao et al. and claimed 
to be unwarranted. Also, the method outlined above is 
a more direct test of the claims of Zhao et al. 

What do we expect from analyses involving such a 
spectral index, if the claims of Zhao et al. are true? First 
of all, we would expect that n^ <10 ° and nf >10 ° converge 
to the values found in runs 1 (or 5) and 2 in the previ- 
ous analysis, respectively, if we do a run over the range 
I = 2 — 220. Further, we would expect that the values 
of R found over such a multipole range would be greater 
than those values found in runs 3 and 6 in the previous 
analysis. Finally, since we are effectively introducing a 
new parameter, we would expect a significantly better fit 
to data - that is, the likelihood values for the ML model 
of this analysis should be well above the likelihood val- 
ues of the ML models of runs 3 and 6 in the previous 
analysis. 

We will vary the same parameters as in the previous 
analysis, both with and without r. In addition, we will 
do runs where we fix R (to see if this makes 100 and 
n s >10 ° converge better to their expected values), and 
runs where we fix n e s <10 ° and 7ig >10 ° (to see if R con- 
verges better to its expected value). When fixing R and 
rig 100 , we fix them to their full-dimensional ML value 
found in run 1, and when fixing rij >10 °, we fix it to its 
full-dimensional ML value found in run 2. 

We will also do a full analysis, similar to the one in 
section HVl the only difference being the use of the step- 
like spectral index. We only do this analysis for the 7- year 
data, since the analysis in section ITVl was only carried out 
for this data set. 

The data used are again the WMAP 5-year and 7-year 
TT,TE,EE, and BB data, and again we do all runs 
twice, once for each set of best-fit parameters. 
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Parameter ML value One-dimensional peaks and MCI limits 

y 68%f, 95%t 
A 68%4, 95?U 

5-year 7-year 5-year 7-year 

Run 1 (1 = 2— 100, varying n s , A s , and R) 
n. 1.038 1.036 1.047 J;™; £™ 1-049 £££ J^g 
log(10 10 A s ) 3.066 3.067 3.047 3.050 ^J;^ 

g 0.104 0.107 0.069 S °;ggg 0-138 g;ggg; g;^g 

Run 2 (i = 101 - 220, varying n s , A a , and H) 
n s 0.925 0.953 0.961 gg3 ^ 828 0.979 914 | .846 
log(10 10 A 3 ) 3.291 3.234 3.203 ^ 3.489 3 m 3.316, 3.455 

/? 028 245 n 00c 1 044, 2.050 n 1 1 , 1.144, 2.161 
u.uzo u.z^to U.OJJ 0.000, 0.000 »- ilz 0.000, 0.000 

Run 3 (I = 2 - 220, varying n s , A s , and 7?) 
^ 0.976 0.991 0.998 Ejjgjliwj 1-010 ^ 
log(10 10 A s ) 3.185 3.158 3.136 ^ 3.116^;^? 

g 0.001 0.021 o.ooo g;jgg; g;gjg o.ooo g;p^; g:ggg 

Run 4 (£ = 2 — 220, varying n s and A s ) 

n. i.ois 1.024 1.(111 ;,;!];::;; 1.021 

log(10 10 A a ) 3.098 3.089 3.098 IZ] jl™ 3.092 g;^g; |;gg* 
Run 5 [£ = 2 — 100, varying r, n s , v4 s , and i?) 
t 0.0975 0.0941 0.1021 g;^ 7 ; g;^| .0937 g;^?; ™ 
1.076 1.063 1.095 ^ 9 7 1.084 ^**;^ 
log(10 10 A s ) 3.025 3.041 2.982 f™- 3.028 ^ 

E> n if? n i /ic n ono 0.341, 0.520 n -, oc 0.293, 0.447 

R 0.167 0.146 0.20 8 044 ; 000 0.136 038 ; 000 

Run 6 (£ = 2 — 220, varying r, n s , v4 s , and R) 
r 0.0876 0.0876 0.0869 ™ ™ 0.0879 Sio?!?: S 
0.980 0.993 1.008 J;"**; J;!™ 1.020 J;°|«; J;™ 
log(10 lf %) 3.191 3.166 3.139 g;™; ^ 3.134 ^f 7 ; 

g 0.002 0.022 0.000 g; gg; g;ggg 0.015 g;^g; g;ggg 

Run 7 = 2 - 220, varying r, n«, and A s ) 
r 0.0882 0.0912 0.0942 <™ g;J™ 0.0927 g;^|; g;»« 

1.020 i.o28 i.o2i i:SoI: S:iIS i.o3o S:2£ 

log(10 10 A s ) 3.108 3.101 3.101 ^;^ 3.094 

TABLE III: The unmarginalised ML points and one-dimensional peaks and MCI limits resulting from a likelihood analysis 
using both 5-year and 7-year WMAP data for various £ ranges, and with various parameters. The best-fit values for the fixed 
parameters are the full-dimensional ML values taken from table [TJ MCI limits are given with both 68% and 95% confidence. 



1. Results 

The results of our analyses are shown in tables |V] (our 
best- fit values) and IVII (the best-fit values of Zhao et al.). 
The results of the full analysis are shown in table I VIII 
Finally, we show the comparison of likelihood values for 
the ML models in tables I Villi (our best-fit values) and 
HXl (best-fit values used by Zhao et al.). 



2. Discussion 

From tables [V] and IVII we see that at least two of our 
expectations have not been met, or only partially: We 



find no higher values for R when using a step-like spec- 
tral index than when we assume the spectral index to 
be constant, and under no circumstances does the value 
of R reach the value it gets in runs 1 and 5 from the 
previous analysis, or anything near it. Further, even 
though n. <1Q0 and n £ s >100 behave somewhat as expected, 
in that n e s <10 ° > nl >wo (so that the spectrum is bluer for 
I < 100 than for I > 100), their distributions lie closer 
than the distributions for n s from runs 1 and 2 (compare 
run 8 with runs 1 and 2). 

Comparing the likelihood values of the ML models, we 
see from tables I Villi and IIXI that our third expectation 
has not been met either: Including a step-like spectral 
index only marginally raises the maximum likelihood, 
which suggests that allowing for an ^-dependency in the 
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Parameter ML value One-dimensional peaks and MCI limits 

y 68%f, 95%t 

5-year 7-year 5-year 7-year 

Run 1 (1 = 2— 100, varying n s , A s , and R) 
n. 1.042 1.046 1.063 ^?;^? 1-062 S;Jgg; J;^* 
log(10 10 A s ) 3.042 3.036 2.997 g;* 7 g; g;?g* 2.974 *;™;°;™ 

g 0.120 0.129 0.140 giggg; giggg 0-132 g;g||; g;gg 

Run 2 (i = 101 - 220, varying n s , A a , and H) 
n. 0.921 0.960 0.965 J;°g£ J;°g* 0.994 J;gg?; ^| 
log(10 10 A 3 ) 3.282 3.207 3.210 |;g£J;«* 3.128 *;2"; ^ 
g 0.046 0.365 0.000 gjggg; jjggg 0-582 g;jgg : g;ggg 

Run 3 (£ = 2 - 220, varying n s , A s , and 7?) 
n. 0.976 0.997 0.998 J;ggg; J;g|g 1-012 J;°£; £;°jjg 
log(10 10 A s ) 3.169 3.132 3.124 *;£?; £j™ 3.098 

g 0.009 0.043 0.000 g;5ig; g;ggg Q.Q38 g;jgg : g;ggg 

Run 4 (i = 2 — 220, varying n s and A s ) 

7, 1 n91 i n oi i n99 1.040, 1.059 i r.01 1.053, 1.072 

n s 1.UZ1 l.UOl l.UZZ 1.001, 0.982 J-.UOX i.on, 0.991 

log(10 10 A a ) 3.075 3.062 3.073 g;^; 3.057 |ffjjjjg 

Run 5 (£ = 2 — 100, varying r, n S) A s , and i?) 
t 0.0965 0.0966 0.1002 g^g; g;J«g 0.0960 g^gg; g;Jgg? 
1.069 1.060 1.081 };J£; l f 6 l 1.078 
log(10 10 A s ) 3.017 3.027 2.977 g;^; 3 2 f 9 l 3.006 *;£»; g;?g* 

D n i co n i co n am 0.344, 0.509 n i i c 0.285, 0.445 

R 0.158 0.152 0.20 1 059 ; 000 0.116 032 ; 000 

Run 6 (£ = 2 — 220, varying r, n S) j4 s , and R) 
r 0.0869 0.0886 0.0924 ™j g; g 9 7 g 0.0912 g;Jggg 
0.976 0.997 1.005 g,;g^ : ;g^ 1.01lS;S 
log(10 lf %) 3.175 3.137 3.128 g;ggg; g;g« 3.117 g;^ ; g;gi? 
R 0.004 0.041 0.000 g;ggg; g;ggg 0.040 g;^g; g;ggg 

Run 7 (I = 2 - 220, varying r, n«, and A s ) 
r 0.0912 0.0932 0.0905 g; ?gg; g; ^ 0.0917 ™ g; gg 6 7 
m 1 nos i n-^ 1 n97 1 046 > 1 067 i oon 1 057 > 1 078 

7l s 1.UZ.5 l.UOO l.UZ/ 1.003, 0.983 J..U3I 1.013, 0.992 

log(10 10 A s ) 3.085 3.067 3.086 g;gg; g;gg 8 3.072 g;^; 

TABLE IV: The unmarginalised ML points and one-dimensional peaks and MCI limits resulting from a likelihood analysis 
using both 5-year and 7-year WMAP data for various £ ranges, and with various parameters. The best-fit values for the fixed 
parameters are the same as those used by Zhao et al. MCI limits are given with both 68% and 95% confidence. 



spectral index serves little purpose, as it means introduc- 
ing a new parameter which does nothing to improve the 
fit of our model to the data. 



C. Data analysis conclusions 

Our analysis seems to suggest that the claims concern- 
ing an ^-dependent spectral index made in [5[ and [6| does 
not hold up under closer scrutiny. It seems there is nei- 
ther any reason to abandon the assumption of a constant 
spectral index, since the distributions do overlap for two 
different £ ranges and since splitting up the spectral in- 
dex as above gives no better fit to data, nor that this 
assumption leads to a negligence of signs of gravitational 



waves in the CMB. 

We saw that in some of the 2 < I < 100 runs, we found 
a lower limit on R. However, we noted that this lower 
limit was almost zero, and that it was just a ltr lower 
limit anyway. Our analysis gave lower R values than the 
analysis in the papers by Zhao et al., which probably is 
mostly due to our use of the exact likelihood functions 
from WMAP in contrast to the approximations made in 
those papers. 

As we mentioned in our comments on the papers by 
Zhao et al., we do not agree that searches for gravita- 
tional waves should be limited to a smaller set of data. 
As long as our cosmological parametrisation is reason- 
able, the Bayesian framework will automatically ensure 
that any signs of gravitational waves will be detected. 
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Parameter ML value One-dimensional peaks and MCI limits 

y 68%f, 95%t 
A 68%4., 95%4 

5-year 7-year 5-year 7-year 

Run 8 (varying nf <1()0 , nf >10 °, A a , and R) 

iwr" l.ois ;°l 7 ; i-02i 5:gg; 

nf 100 0.983 0.987 1.004 ;° 3 7 7 ; J;gjg 1.016 
log(10 10 A) 3.172 3.165 3.124 l^; 3 ; 2 7 4 3.117 I;™; 3 ; 22 ? 

J? 0.004 0.004 o.ooo g;jgg; g;gg; o.ois g;jgg; g:ggg 

Run 9 (varying r, nf <10 °, n e > 100 , A s , and g) 

r 0.0824 0.0926 0.0874 ™ 0.1286 q Q900 ™ ™ 

£<100 1 004 1 q 27 1 Q21 1.093, 1.160 -, Q44 1.111, 1.181 

n s l.UWi 1.UZ1 .956, 0.888 J--U*a 0.962, 0.896 

„£>ioo n qoi n QQ7 1 rw 1 048 ' 1 090 1 rwn 1 063 ' 1109 
n 3 U.981 0.997 1.UU4 . 978i .950 L - 0A0 0.989, 0.957 

log(10 10 A s ) 3.176 3.169 3.138 ^ %™ 3.124^;^ 

R 0.003 0.017 0.000 g;^; gjggg 0-000 qogp; qgoo 

Run 10 (varying r, nf <10 °, nf >10 °, and As) 

r 0.0894 0.0924 0.0906 ™ ™ 0.0929 ™j ™ 

£<ioo 1 n2C - 1 n oq 1 n31 1 096 ' 1 162 1 025 1 105 ' 1178 

"s l.UZO X.UO» X.UOX 0.960, 0.893 X.UiJ 0.966, 0.898 

£>100 , n22 1 o 2 j q 21 1.049, 1.073 m4 1.055, 1.080 

"s l.UZZ X.IX4X 0.999, 0.975 X.UOI 1.005, 0.980 

log(10 10 ^) 3.105 3.093 3.109 f^tlM 3.106 jZ 
Run 11 (varying r, A s , and i?) 
r 0.0726 0.0830 0.0722 « ™ 0.0825 g;gg«; ™ 

log(10 10 A s ) 3.257 3.230 3.260 3 ; 3 J 2 3.226 3 : 27 ? 

0.011, 0.029 n nnn 0.018, 0.043 
0.000, 0.000 U-UUU 0.000, 0.000 

Run 12 (varying ni <10 ° , n l > 100 , and A a ) 

!<<ioo Yom 1 022 1 013 1 ' 076, 1 141 1 019 1 ' 080, 1 153 
n a i.uu» l.uaa x.uxo 0.943, 0.882 x.ux» o.952, 0.894 

£>100 lnl7 j q 24 j ni o 1-038, 1.060 1 n22 1.046, 1.068 

n s l.UK 1.UZ4 I.UIO 0.994, 0.973 l.UZZ 1000> . 979 

log(10 10 A s ) 3.101 3.089 3.100 



R 4 x 10- 5 2 x 10- 4 0.000 g'LV n mn 0.000 °' 018 ' 043 



3.148, 3.189 o n oo 3.137, 3.179 

3.058, 3.012 O.U58 3,045 3 000 



Run 13 (varying A s and R) 


log(10 1( M s ) 
R 


3.271 

8 x 10" 


nor o 1 ?/! 3.280, 3.288 
O.Ai^-k g 262 3 254 

6 2 x 10~ 5 000 012 ' 031 
Z X 1U u.uuu 0.000, 0.000 


3.223 
0.000 


3.232, 
3.215, 
0.019, 
0.000. 


3.240 
3.207 
0.044 
0.000 



TABLE V: Unmarginalised ML points and one-dimensional peaks and MCI limits resulting from a likelihood analysis with an 
^-dependent n s , using the best-fit values derived in this paper. Results from the 5-year and 7- year analyses are shown, and all 
runs were over the range I — 2 — 220. MCI limits are given with both 68% and 95% confidence. 



The analysis performed in this paper brings us one step 
closer to saying that our parametrisation indeed is rea- 
sonable. 

In any case, we have seen that the only indications of 
deviations from the standard analyses are very weak, only 
on the level of one confidence interval, and any claims 
made based on such weak detections become little more 
than speculation. As mentioned in the papers by Zhao et 
al., ongoing experiments such as Planck [18| and QUIET 
19], among others, should make us much better suited to 
answer such questions as those addressed in this paper. 



VI. THEORETICAL CONSIDERATIONS 

It seems plausible that at least some of the motivation 
for the papers [f| and [(| comes from the work of Gr- 



ishchuk on relic gravitational waves. He discovered that 
in any accelerating universe, existing gravitational waves 
will be amplified as long as their wavelengths are long 
enough (20[. It is also notable that Grishchuk, since the 
mid-90s, has been a voiced opponent of the inflationary 
paradigm, which also may have been a motivation for the 
data analyses described above. 

In Q , density perturbations are evolved through three 
stages: an initial, possibly inflationary stage, a radiation 
dominated stage, and a matter dominated stage. The 
amplitude of the density perturbations are then com- 
pared to the amplitude of the gravitational waves in the 
matter dominated stage (which were derived in an ear- 
lier paper, see [Hj]), an d Grishchuk finds that the density 
perturbations do not dominate; instead, the gravitational 
waves are the dominant contribution to the CMB. This 
is contrary to the predictions of the standard derivation, 
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Parameter ML value One-dimensional peaks and MCI limits 

y 68%t, 95%t 
68%4, 95%4. 

5-ycar 7-year 5-year 7-year 

Run 8 (varying nj <10 ° , n e > 100 , A s , and R) 
~^^> 1^020 1*14 1.011 ;°t Hfr 1-032 ^ 
nf ™ 0.985 0.993 1.006 ;°% ;°™ 1.012 ;°«; 
log(10 10 A s ) 3.151 3.140 3.120 ^J; 22 9 7 3.103 !;J£; |;™ 

g 0.006 0.019 o.ooo g;ogg; g;ggg o.ois g;£gg; gjjg 

Run 9 (varying r, nf <10 °, nf >10 °, y4 s , and R) 

r 0.0888 0.0904 0.0905 j™ ™ 0.0920 ™ ™ 

«<100 i gig i 02fi i q 35 1.100, 1.167 1 041 1.111, 1.177 

' l s I.UIO l.UZO l.UOO o.960, 0.891 l.UfH .970, 0.900 

£>ioo 984 Q 996 i nn 5 1 046, 1.092 1 m4 1.059, 1.103 

' l s u.trai u.yyo l.UUO . 9 76, 0.948 1 - UM 0.985, 0.956 

log(10 10 A s ) 3.163 3.141 3.121 1;™S 3.105 £££ 

0.002 0.020 0.000 g;g 7 0.015 g;^g; °;ggg 

Run 10 (varying r, nf <10 °, 7if >10 °, and A s ) 

0.0927 0.0907 0.0949 g;Jjj£ 0.0938 °; ?^; g;J™ 

£<ioo 1 n o 2 1 042 1 040 1 ' 104, 1 177 1 043 1 n1, 1184 

n s 1.U.5Z 1.U1Z 1.U1U 0.965, 0.897 Lu * i 0.971, 0.907 

<->100 j n2 t- 1 Q 4 1 noi 1 054, 1.081 j ^ 1.062, 1.089 

n s l.UZO l.UOl l.UOl 1001j 0.978 10 11, 0.986 

log(10 10 A s ) 3.084 3.062 3.076 |^ 3.061 ^ 4 

Run 11 (varying r, j4 s , and R) 

aTvtTc n noe-7 nnin 0.0885, 0.1023 nnM . 0.0973, 0.1101 

r 0.0745 0.0867 0.0711 a059g ; . 04 59 °- 0824 0.0708! o.osso 

10^ ^ q OAR Qine Q OAA 3.272, 3.301 Q 1 no 3.218, 3.244 



log(10 iu A s ) 3.246 3.198 3.244 



R 1 x 10- 4 5 x 10- 5 0.000 g'g, 1 ,, 1 / g'gg 8 , 0.000 °' 021 ' 050 



3.214, 3.186 
0.011, 0.028 

O.OOO, 0.000 0.01)0, 0.000 



Run 12 (varying n l < lm , < >10 °, and .4 S ) 



0.953 



< <iuu 1.023 1.049 1.001 
n e > 100 1.022 1.036 1.0 1 8 „.,,„„ 
log(10 10 A) 3.073 3.052 3.070 l™> ^ 3.162^;^° 



1.153 n Q9Q 0.980, 1.028 

0.889 U.yzy o.879, 0.827 

1.065 n Q77 0.984, 0.992 

0.973 U - y ' ' 0.967, 0.958 

3.120, 3.171 o 1eo 3.186, 3.208 



Run 13 (varying A s and 7?) 



log(10 10 A s ) 3.263 3.199 3.263 3 3 ^XZ 3.197 ^°g; ^ 
R 4 x 10~ 5 4 x 10~ 5 000 °' 011 ' 027 000 °' 020 ' 048 

£^ 2 X 1U _ X 1U U-UUU 0.000, 0.000 U-UUU 0.000, 0.000 

TABLE VI: Unmarginalised ML points and one-dimensional peaks and MCI limits resulting from a likelihood analysis with an 
^-dependent n a , using the best-fit values used by Zhao et al. Results from the 5- year and 7-year analyses are shown, and all 
runs were over the range I = 2 — 220. MCI limits are given with both 68% and 95% confidence. 



Parameter ML value One-dimensional peaks and MCI limits 

y 68%t, 95%t 

68%j, 95%4 



Q b h 2 


0.0227 


n f\0-}0 0.0241, 0.0250 
U.U/OO 0.0225, 0.0218 


^DMh 2 


0.1096 


n mc;-? 0.1119, 0.1184 

U.1UOO 0.0987, 0.0920 


e 


1.040 


1 041 1 044 ' 1047 

LO^ 1 1.038, 1.035 


T 


0.0910 


0928 °' 1086 ' °- 1275 

U.uy^O 0.0753, 0.0606 


n e < 100 


1.010 


1 n1 -i 1.097, 1.168 
l' uii 0.949, 0.880 


ni> 100 


0.974 


n 088 1 ' 014 > 1043 

U.yOO o ggg^ 0.947 


log(10 ll, A s ) 


3.166 


0111 3.182, 3.246 
> 5 - 111 3.032, 2.941 


R 


0.011 


009 °' 081 ' 0165 
U.UV.J 0.000, 0.000 



TABLE VII: ML points and marginalised one-dimensional 
peaks and MCI limits for a likelihood analysis for eight param- 
eters (including a step-like spectral index) using the WMAP 
7-year data. MCI limits are given with both 68% and 95% 
confidence. 



Parameters varied 


log£ (5-year) log£ (7-year) 


n s , A a , and R (original analysis) 


284.938 


134.985 


n s and A s (original analysis) 


284.244 


134.622 


ni <100 ,n l > 100 , A 3 ,andR 


284.972 


135.044 


nf <10 °, ni> 100 , and A s 


284.243 


134.622 


A s and R 


275.896 


131.768 


t, n s , A s , and R (original analysis) 


284.981 


135.115 


r, n s , and A s (original analysis) 


284.386 


134.832 


r, n e s < 100 , ni >wo , As, and R 


285.027 


135.166 


r, n e s< 100 , ni >wo , and A s 


284.392 


134.816 


t, As, and R 


276.034 


131.739 



TABLE VIII: Log(likelihood) values of the ML points of var- 
ious runs from both the first and second analyses using the 
best-fit values obtained in this paper. The Grange is always 
2-220. 
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Parameters varied 



n a , A 3 , and R (original analysis) 


285.083 


134.921 


Ti s and ^4 S (original analysis) 


284.385 


134.517 


n e s <wo , n e s >lm , A a , and R 


285.205 


135.033 


nf <10 °, nf 100 , and A s 


284.384 


134.490 


A s and i? 


275.849 


132.922 


r, n s , A s , and i? (original analysis) 


284.113 


134.929 


t, n s , and A a (original analysis) 


284.482 


134.565 


t, n{< 100 , n{ >wo , As, and R 


285.274 


135.090 


r, n e a <100 , n e > 100 , and A a 


284.477 


134.566 


t, As, and R 


276.077 


132.923 



log£ (5-year) log£ (7- year) 



TABLE IX: Log(likelihood) values of the ML points of various 
runs from both the first and second analyses, using the best-fit 
values used by Zhao et al. The Grange is always 2-220. 



in which density perturbations dominate. 

In addition to the above papers, Grishchuk authored 
[22| and (23|, in which he presented other arguments 
against the inflationary paradigm. 

The above papers were criticized in [241, to which Gr- 
ishchuk responded with 



251 ] . The argument was ended 



with [26[, though it seems to us to be a bit unclear as 
to what the final conclusions were. After these papers, 
Grishchuk turned to the question of initial conditions in 
[27j , the derivation in which again gave a different result 
than the standard one. This treatment was criticized in 

M- 

Overall, there seems to be little agreement on who, 
exactly, are correct in their treatment. We will therefore 
go through some points of disagreement and try to offer 
our own thoughts on them. 



A. The continuity of n 

Working in the synchronous gauge with the metric 
-a 2 , 9oi = 0, gij = a 2 ((1 + HQS, 



ffoo 



hin 



(6) 



where n is the wavenumber, h is the scalar perturbation, 
and hi is the longitudinal-longitudinal perturbation (so 
that Q becomes the amplitude of the perturbations) , Gr- 
ishchuk aims to find solutions for h and hi, along with 
general variables that appear in the Einstein equations, 
for all the three stages described above. For these stages, 
the scale factor may be written as 



1+/3 



(7) 



where the value of j3 depends on which stage we are con- 
sidering, and r\ is conformal time. 

In addition, Grishchuk defines certain auxilliary vari- 
ables: 



a 



1 



(8) 



where overdots signifies derivatives with respect to con- 
formal time. 7 then becomes a 'smallness parameter' 
during inflation. 

He then goes on to solve the Einstein equations for all 
variables in each of the three stages. It turns out that 
all variables can be expressed in terms of h and hi for a 
given stage. It is therefore sufficient to find solutions for 
these variables. 

At the initial (possibly inflationary) stage driven by a 
scalar field, so that cf>o is the homogeneous part and the 
total field is written as 



iQ, 



(9) 



all variables are expressible through h, and so, Grishchuk 
needs only find a solution for this variable. To this end, 
he introduces the variable [i, defined through 



(10) 



where u = h + orfh. During the initial phase, [i becomes 

H= (m/) 1/2 (AiJ j8+4 (n»;)+AaJ_ os+ i ) (rM 7 )) (11) 

where J denotes the spherical Bessel functions, j3 is the 
variable introduced in eq. (J7]), and A\ and A-i are con- 
stants to be determined by the initial conditions . It can 
be shown that h can then be written as 



a J no a 



(12) 



where Ct is a constant that arises because of the degrees 
of freedom that remain in the synchronous gauge, and 
will be fixed by a choice of coordinate system. 

Grishchuk finds similar solutions for h (and h{] during 
the radiation-dominated and matter-dominated stages, 
with three new constants arising at each stage. Finally, 
he aims to join these solutions, claiming that h, h, and 
hi are continuous over the transitions, and that in addi- 
tion, fi is continuous over the transition from the initial 
to the radiation dominated stage. Using this, he joins 
the three stages, and by choosing the coordinate system 
that moves with the matter fluid, he manages to express 
the unknown constants at the matter dominated stage 
in terms of the constants A\ and A2, which are to be 
determined by initial conditions. He further makes some 
approximations which reduce /1 from eq. (|11[) to 



Km 



2^+*r(/3 + f) 



(nr?i) 



/8+1 



(13) 



where T denotes the gamma function, and 7/1 is the time 
of the transition between the initial stage and the ra- 
diation dominated stage. This means that all unknown 
constants can be related to A\. 

Using this solution, Grishchuk works out the observ- 
able consequences of it, and comparing them to the cor- 
responding results for gravitational waves, he finds that 
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for the CMB quadrupole, the gravitational wave contri- 
bution should actually dominate over the contribution 
from the density perturbations. This does not fit well 
with the standard picture, in which density perturbations 
should be the dominant contributor. 

In [24j . it is pointed out that the variable p can not 
be continuous across this transition, but rather, that the 
combination — ^= is continuous. The authors then claim 
that assuming \i to be continuous is what leads Grishchuk 
to deviate from the standard result. 



Then, in [25j, Grishchuk agrees that /i indeed is dis- 



continuous over the transition, but that this was a typo 
in the original paper, and that he never actually used this 
assumption in his derivation. Finally, in [26j , the authors 
uphold their claim that this assumption was used. 

Grishchuk never responded to this, and so it is hard to 
tell whether he finally agreed with Martin & Schwarz or 
not. One may, however, verify who is right, by looking 
at eq. (|T3|) (as also pointed out by Martin and Schwarz). 
Here, Grishchuk approximates ji at the transition be- 
tween the initial and radiation stages. However, as we 
have already mentioned, /i is discontinuous at the tran- 
sition, which means that one must choose whether one 
wants to use the value of fx at the beginning of the transi- 
tion, or at the end of it. Grishchuk has earlier set 7 = 2, 
which means that he intends to be at the end of the tran- 
sition when he is doing the joining. However, eq. (|13l) is 
an approximation of eq. (jlll) . which is valid during the 
initial phase only. Thus, eq. (ITcH) is the value of \l at 
the beginning of the transition, not the end, as it should 
be. This means that Grishchuk did implicitly use the as- 
sumption of a continuous as otherwise, eq. (1131) would 
not be valid (since we are at the end of the transition). 

Since we have 



/i 



(m + <t) 



(14) 



where a is some arbitrarily small time interval so that 
r/i — a is the beginning of the transition while r\\ + a 
is the end of the transition, we find that the value of fi 
at the end of the transition, which is what we are really 
after in eq. (TT5T) . should be 



Km + a) ~ /i(r/i - a) 



(15) 



where 2 is the value of 7 at the radiation stage, while 7, is 
its value at the initial stage (typically very small). Thus, 

eq. (fl3|) should really be multiplied by J — - This intro- 



duces an amplifying factor which propagates through the 
rest of Grishchuk's derivation, again yielding the stan- 
dard result. 

Thus, it seems that at this point, Grishchuk's treat- 
ment does not hold up under closer scrutiny. 



B. Initial conditions 

Martin & Schwarz adressed some of the other com- 
plaints Grishchuk voiced about the inflationary theory, 
namely, the usage of the constancy of the curvature per- 
turbation C (Grishchuk claims that this variable is equal 
to zero, and thus it is meaningless to utilize its constancy 
in theoretical derivations) and the large amplification of 
the density perturbation ^ from the inflationary stage to 
the radiation dominated stage. Both these issues seems 
to us to be cleared up: Grishchuk seems to agree in 25 1 
that C is both constant and nonzero, and, though it seems 
he disagrees with the exact word use, he also agrees that 
there is an amplification of from the inflationary stage 
to the radiation stage. However, even though ^ is am- 
plified, it does not mean that it necessarily dominates 
over gravitational waves - in order to say anything about 
which component dominates, we must look at the initial 
conditions. 

This is exactly what Grishchuk does in [27{ ■ In this 
paper, he goes through the initial conditions of the early 
Universe both for gravitational waves and density pertur- 
bations, and claims that the standard way of quantizing 
the curvature perturbation £ is incorrect, since (as he 
shows) the expectation value of the normalized pertur- 
bation, called (, is not h/2 as it should be for a proper 
harmonic oscillator; it is proportional to the 'largeness 
factor' 1/7,. This, Grishchuk says, is due to the fact 
that the quantum states that are used for quantization 
in the standard derivation are not ground states; they are 
so-called 'squeezed states'. By doing a Bogliubov trans- 
formation of these states, he ends up with new states 
that satisfies the ground state condition. With these ini- 
tial conditions, he finds that the spectrum of £ at horizon 
crossing is comparable to, not hugely amplified compared 
to, the spectrum of gravitational waves. 

His treatment was criticized by Lukash in (29j j , in which 
the author says that Grishchuk makes an error in the 
normalization of £, and it therefore becomes an error 
to demand the expectation value of this variable to be 
that of the ground state of an harmonic oscillator. In 
short, Grishchuk blamed the states, when he should have 
blamed the normalization of 

Lukash makes some other points as well to back up his 
criticism, but this will do for our purposes. We do agree 
with Lukash that the normalization of £ in Grishchuk's 
paper seems inconsistent with what he himself says: £ 
is normalized in exactly the same way as h (the tensor 
perturbation), namely: 



c 



pi 



(16) 



where Mpi is the Planck mass, and do is the value of the 
scale factor at some initial time at which to fix the initial 
conditions (with 70 being the corresponding value of 7,. 
This value changes very little during inflation, so we have 
70 w 7j). However, he also says that when moving from 
the treatment of gravitational waves to density perturba- 
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tions, one should always do the replacement a — ¥ a^yl. 
This means that the normalized curvature perturbation 
really should contain an extra factor -y/To, which gets rid 
of the factor when calculating the expectation value 
of the £ spectrum. Therefore it seems to us that Gr- 
ishchuk's claims about the initial conditions have been 
rebutted by Lukash. 

VII. CONCLUSIONS 

In this paper, we have studied various claims concern- 
ing primordial gravitational waves: We first looked at 
the claims made in [f| and @ concerning the detection 
of a nonzero contribution of gravitational waves in the 
WMAP 5-year and 7-year data, by limiting the multipole 
range to I = 100. After first noting that limiting the mul- 
tipoles thus is neither more correct (as claimed by Zhao et 
al.) nor optimal in these data analyses, we repeated their 
analyses. We found that by using CosmoMC and the of- 
ficial noise values instead of approximated values for the 
noise, the gravitational wave contribution lessened, and if 
lower limits for R were found, they were typically small. 
We also found no reason to believe their claim that it 
is an error to assume a constant spectral index for all 
multipoles, as the distributions for n 3 overlapped for the 
t = 2 - 100 and I = 101 - 220 analyses. 

Nevertheless, we tested this claim further by imple- 
menting an ^-dependent spectral index, allowing n s to 
have two different values above and below £ = 100. Our 
findings were not what we expected them to be if the 



claims by Zhao et al. are correct: The values of n s <wo 
and rig 100 did not converge to their expected values, 
and even with such a step-like spectral index, the gravi- 
tational wave contribution did not rise to any significant 
level. 

We did all of the above analyses twice; once where 
the parameters that were not varied were given the same 
values as by Zhao et al., which were the mean values 
from the WMAP analyses, and once where we used val- 
ues from our own initial analysis where we allowed for 
gravitational waves but used the relation n t — n s — 1 in 
accord with the statements of Zhao et al. For both sets 
of analyses, the results were qualitatively the same, so 
the conclusions seem to be independent of which best-fit 
values one prefers to use. 

Finally, we went through some of Grishchuk's ear- 
lier claims concerning primordial gravitational waves, as 
these were conceived to be at least part of the motivation 
behind the paper by Zhao et al. We reviewed the contro- 
versy that followed Grishchuk's claims, and offered our 
own thoughts on the matter. In short, it seemed to us as 
Grishchuk's claims were firmly rebutted and that he has 
not responded adequately to these rebuttals. 

The search for gravitational waves in the CMB will 
undoubtedly continue, and satelites such as Planck will 
surely shed further light on this issue. Just as surely, 
there are bound to be surprises in the next generation 
of CMB data. However, based on the data we have to- 
day, the gravitational wave contribution to the CMB fits 
nicely with the Standard Model and the inflationary the- 
ory. 
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